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Three subjects are considered here: the relativistic hydrodynamics equations for a boost-invariant 
expanding fluid; the fuzzy bag model for the pressure which recently appeared in QCD phenomenol- 
ogy; and the early space-time evolution of the QCD matter, drawn from model studies, which can 
also be expected to arise in realistic fluid dynamics relevant to heavy ion collisions at LHC. 
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I. INTRODUCTION 

The idea of exploiting the laws of hydrodynamics to describe the expansion of the strongly interacting matter 
created in high energy collisions was pioneered by Landau in 1953 [1 . In recent years, data from the Relativistic 
Heavy Ion Collider (RHIC) provided striking evidence for a strong collective expansion that is in good quantitative 
agreement with hydrodynamic predictions.^ Moreover, these data emphasized the exciting possibility that the QCD 
matter created in the collisions shows properties of a near-perfect fluid. Once well-calibrated ideal fluid dynamical 
benchmarks have been established under appropriate experimental conditions, deviations from perfect fluid behavior 
can be used to explore transport properties, such as viscosities and heat conductivities of the QCD matter. If 
unambiguously extracted from data, they are prime candidates for the next compilation of the Particle Data Group 
[3]. Such efforts are among those that define the present forefront of research in heavy- ion collision dynamics. 

The results at RHIC also indicate that the matter created in heavy ion collisions behaves like a strongly coupled 
liquid rather than a weakly interacting plasma of quarks and gluons. In this case, the MIT bag model that effectively 
includes strong interactions via a bag constant can be of use as an equation of state. Recently, the bag model has 
taken an interesting turn with the detailed analysis [4j of the lattice data [5^ . The surprise of this analysis is that at 
zero chemical potential it reveals a term quadratic in temperature as the leading correction to the ideal gas term in 
the pressure.^ The recent data [7H9] also indicate that this is generic: that the equation of state of the original bag 
model gets modified as ^ 

p{T) = a{T* - aT'^) - B , with r„,„ < T < r„„, . (1) 

As noted, the novelty is the term, while the remaining two are the standard bag terms, with B a bag constant 
and a a parameter. A common choice is to take a from perturbation theory up to one loop order. r„,i„ is close to a 
critical temperature Tc (or some approximate Tc for a crossover). A small difference between Tc and T,„i„ may vary 
with the model. T^^^ is set by perturbation theory such that to leading orders it is applicable only for temperatures 
higher than T^^^. In fact, there are arguments in the literature [10] that in QCD long perturbative series (or the UV 
renormalons) result in the so-called quadratic corrections. From this point of view, the term in ([l]) is nothing 
else but an example of the quadratic correction. In the dual (string) description, these quadratic corrections are 
parameterized in terms of the metric in extra dimensions and do result in the term for the pressure jllj . 

What is the physical significance of such a modification? In the present paper, we look into how it affects the early 
space-time evolution of the QCD matter and why it may be relevant to some attempts to extract dissipative transport 
coefficients from data. 

In general, the 3 -t- 1-dimensional hydrodynamic equations are very complicated and require powerful numerical 
simulations. Here, we consider a simplified model which, however, is a good approximation to the early time evolution 
of the QCD matter. It incorporates several characteristic features of heavy ion collisions, while lacking many of the 
technical complications of the most general case. Moreover, it allows us to find analytic solutions. This is the fluid 
dynamical model proposed by Bjorken [12]. 

In this 1 -|- 1-dimensional model, introducing the Milne variables proper time t = \Jt'^ — and space-time rapidity 
^ = arctanh(z/t), boost invariance simply translates to requiring all hydrodynamic variables to be independent of 



^ For reviews, see 0. 

^ Phenomenologically, such a term was first suggested in [6]. 
^ Rob Pisarski called it a fuzzy bag model for the pressure [4]. 
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rapidity. As a result, the hydrodynamic equations become exceptionally simple [2- To first order in the gradient 
expansion, one has 
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given by their Navier-Stokes values. Here 77 and ( are the shear and bulk viscosity coefficients, respectively. 

In order for the equations ^ and ([3| to be closed, one has to specify an equation of state p{e) and expressions 
relating 77 and C to e. The simplest case to be considered is that of a perfect fluid, 77 = C = 0, with a bag model type 
equation of state, eq.Q at a = 0. If so, then the evolution equation ^ can be solved analytically. One gets 



(4) 



with the initial condition T(to) = Tq. This is the celebrated Bjorken solution 12J which defines the baseline on top 
of which dissipative and fuzzy bag effects have to be established. Certainly, deviations must be sufficiently small for 
the gradient expansion underlying dissipative hydrodynamics to be valid. 



II. PERFECT FLUID 



In the absence of dissipative effects, the entropy per unit of rapidity is a constant of the motion [T^]. Therefore the 
entropy density is 



s(r) = So 
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For the fuzzy bag ([!]), we can still obtain an explicit formula for T{t) by solving a cubic equation. So we get 
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where c = toTo{T^ - i")- The solution Q emerges in the limit a — >■ 0, as expected. 

In Figure 1 we show the proper time evolution of the temperature and energy density for the bag and fuzzy bag 
models. Following our discussion on the fuzzy bag, we choose the initial temperature to correspond to T^^^ and stop 
the evolution at T = T„i„. It is expected on general grounds that r„i„ is somewhere between 0.2 GeV and 0.3 GeV, 
while T,-„„^ is around 1 GeV. This is consistent with the recent lattice data for QCD with almost physical quark masses 
in the interval 0.3 GeV <T< 0.8 GeV [7 . So, we take the initial temperature To to be 0.8 GeV. It is worth noting 
that this value is expected to be reached in the coming years as the LHC goes into its heavy ion program. Given To, 
the initial time can be estimated by using the uncertainty principle |13j : To{E)q ~ 1, where {E)q ~ 3To. This results 
in To = 0.08 fm/c. Thus, we consider the early time evolution where the Bjorken model is a good approximation. The 
value of a is fixed from the slopes for the Regge trajectory of ^(77) mesons and the linear term of the Cornell potential 
[I4] . For Nc = Nf = 3, it is given by a « 0.09 GeV^. Finally, the value of the bag constant is set to (0.2 GeV)'*. 

What we see from this Figure is that the temperature decreases with the proper time more slowly in the presence 
of the term in the equation of state. Note that the bag constant has no effect on the evolution of the temperature. 
Thus, the Bjorken solution is also valid for a conformal fluid. The maximum discrepancy occurring at r = 2 fm/c is 
of order 13%. It rapidly decreases with decreasing t and becomes almost negligible for r < 0.3 fm/c. 

For practical purposes, the expression ^ looks somewhat awkward. A possible way out is to expand it in a series 
and then see where that takes us. If we ignore all higher terms in the small a expansion, which is equivalent to the 
small T expansion, then a final result can be written in a simple form 
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FIG. 1: Left: Evolution of the temperature. Right: Evolution of the energy density. In both the cases, the lower blue and 
upper red curves correspond to the MIT and fuzzy bag models, respectively. The dots come from the approximation ([7|. 



Here the second term represents the leading a-correction to the Bjorken solution. Since its value is positive, the 
slowing down of the temperature decrease occurs (at least for small r). In Figure 1 we have plotted the proper time 
dependence of the temperature as it follows from Q. Remarkably, with our initial conditions, the difference between 
([g]) and Q is negligible in the interval 0.08 fm/c < r < 2fm/c. 

Finally, let us discuss the evolution of the energy density. In the fuzzy bag model it is given by e = SaT* — aaT^+i?.^ 
Because of the minus sign in front of the term, one would expect that the energy density will fall faster in the 
presence of the term. In reality, the effect of the slowing down of the evolution of T dominates and, as a consequence, 
the energy density falls slower in the fuzzy bag model than in the original one. We illustrate this in the right panel 
of Figure 1. 



III. VISCOUS FLUID 



One of the most fascinating developments in recent years has been the conjecture |15| that there may be a funda- 
mental bound from below on the ratio of the shear viscosity 77 to the entropy density s ^ 

47r s ' 

Given the conjectured existence of the bound, it is natural to ask how closely does the QCD matter (fluid) produced 
at RHIC (LHC) approach the bound. While some methods provide consistent support for a conclusion that the 
produced matter is within a factor of 2 — 3 above the conjectured bound [16^, controversy still remains [T7]. We now 
wish to understand whether the modification of the original bag model will affect the conclusion. 

A. Shear Viscosity or Modified Equation of State? 

In the case of the MIT bag model, the solution of the evolution equation (|2| is known analytically for C = and 
77/5 — const [H] 



T(r)=To(^)% 

where k = 77/s. The second term is a viscous correction, 
time more slowly in the presence of the shear viscosity. 

Now, a problem arises. For r < tq/{\ + fl^To)^/^, the 
What it means is that dissipative effects are large and, as 
more.® To proceed, therefore, we need to pick the initial 




It shows that the temperature decreases with the proper 

function ([9| grows with increasing r. This is unphysical. 
a consequence, the first order theory is not applicable any 
time in such a way that tq > to/(| + |^?o)^''^- Together 



^ For completeness, we have included a formula for tij) in the Appendix. 

We use units where h = c = = 1. 
^ This is an example of reheating. For more discussion, see 0] and references therein. 
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with tqTq = 1/3 13J, this results in the following constraint k < 1/4. Since the value of l/(47r) is approximately 0.08, 
it obeys the constraint. Moreover, it is possible to consider larger values of k (up to a factor of 3). Finally, let us note 
that the situation is opposite for the solution (|6]), where the problem of small r doesn't occur. 

From our discussion so far, when the Bjorken solution Q defines a baseline, the slowing down of the evolution of 
the temperature occurs if either the equation of state gets modified or the shear viscosity coefficient is not set to zero. 
A natural question arises: what is dominant? 

We can gain some understanding of this by plotting the curves given by equations ^ and The results are shown 
in Figure 2. We see that the function (|6| decreases with r more slowly. Indeed, at the beginning of the evolution 



T (GeV) 




FIG. 2: Evolution of the temperature. The green curve corresponding to k = l/(47r) lies between those for k = 0.03 and 
K = 0.12. The dots represent the solution Q with a — 0.09 GeV^. As before, the initial conditions are To = 0.8 GeV and 
TO = 0.08 fm/c. 

process it is comparable with the function ([9| at values of k below l/(47r), while at the end it becomes comparable 
with that at greater values exceeding l/(47r). Also note that for k = l/(47r) the functions are almost indistinguishable 
in the interval 1.2fm/c ^ t < 1.7fm/c. 

In fact, it is easy to understand why it is so. The reason for this is contained in the formulas ^ and The 
bmction ^ goes to \/a/2 as t — > oo, while the function ^ goes to 0. We will return to this issue in section VI. 

Thus, what we have learned is that these effects have the same sign and both are significant in the interval of 
primary interest 0.08 < rj/s < 0.24. We consider the early time evolution 0.08 fm/c < r < 2 — 2.5 fm/c, where 
the Bjorken model is a good approximation. In such a situation, our conclusion is likely to be reliable in heavy ion 
collisions as well. If so, it should be taken seriously by those who will try to extract the value of the shear viscosity 
from the LHC data and shed some light on the conjectured bound. 



B. Attempt of Synthesis 



The formula ([6]) is oversimplified for various reasons. For one thing, the shear and bulk viscosity coefficients must 
be included in the evolution equation. A common way to do so is to set r//s to a constant and C to zero [2]. For 
the window 0.3 GeV < T < 0.8 GeV, it gives a good approximation for the evolution. But this level of accuracy is 
sufficient nether for the next compilation of the Particle Data Group nor for shedding the light on the lower bound 




First, we will discuss the shear viscosity. Like in atomic and molecular physics, the ratio rj/s depends on T. In the 
temperature interval of interest it is an increasing function. Although there is strong evidence for this from lattice 
simulations [20] , the limited data don't allow one to suggest any reasonable parametrization. On the other hand, in 
an alternative approach based on the soft wall metric model of AdS/QCD a simple formula was derived in |21j. For 
T > Tc, it is 



(10) 



where A: is a positive constant and a is given by 0. The reason for this is that the model was also used to derive the 
equation of state jTT]. In contrast to [5T], we have not normalized (10) to k = l/(47r) 



The situation seems to be better with the bulk viscosity. According to ^22j, C can be related with the QCD trace 
anomaly, yielding for high temperatures (/s ^ l/T^. This conclusion is attractive because it agrees with what can 
be drawn from lattice simulations [57. In fact, the lattice data are well fitted by [21] 
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Here m is a positive constant and a is introduced for dimensional reasons. 

This state of affairs, together with the evolution equation ([2]), means that we can combine the viscous terms as 
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where 5 — k — 3m/ (4k). 

There is an obvious question. Can the ratio ^/s, which is dependent of T, cancel the T dependence of r]/ si Note 
that this is not merely a problem with the choice of fitting functions. It is a more general phenomenological problem. 
The ratio "q/s is an increasing function of T, while the ratio C/s is a decreasing function of T. Therefore, if one 
combines these two, a cancellation among "T dependencies" can happen. 

We have no real resolution of this problem. All what we can do is to illustrate it with an example. We do so in 
Figure 3. Here we take k = 1/A at k = l/(47r) as it follows from the results of [51] and [TT]. It is worth noting that 
r]/s ~ 0.06 at T = 0.3 GeV that is of course below the conjectured lower bound on rj/s of [H] but above 0.03 of [TT]- 
With the same values of parameters in ( [Tl] ) as those of [24], we have m « 0.024 at a = 0.09 GcV. Then we make an 
estimate of i5 « 0.024. Surprisingly, it is one order of magnitude smaller than k. We see therefore that f is a slowly 
varying function of T: there is a small difference between the conjectured value of k = l/(47r) and v. The maximum 
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T (GeV) 

0.4 0.5 0.6 0.7 0.8 
FIG. 3: Values of v and rj/s versus T. The horizontal line represents the conjectured lower bound on rj/s. 

discrepancy occurring at T = 0.3 GeV is of order 2%. For comparison, note that it becomes of order 25% for rj/s. 

Having determined the viscosity coefficients, we can now consider the evolution equation. Combining ([ij and (12 1, 
we get 



2„ /„, 1 \ a 



2^2 
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We cannot unfortunately solve this equation analytically but wc can do it numerically. Figure 4 shows the time 

T (GeV) 
0.38 
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FIG. 4: Evolution of the temperature. Left: The curve corresponding to 5 = 0.024 lies between those for S = ±5. The dots 
come from the approximation ( |14[ | a,t 5 — 0.024. Right: The solid curve corresponds to 5 — 0.024. The dashed curve represents 
the approximation (16 1. 



dependence of the temperature for our initial conditions To = 0.8 GeV, tq = 0.08 fm/c. We also take a = 0.09 GeV^ 
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and K — I/Att. In fact, to see significant effects on T{t), the absolute value of the parameter S must be increased by 
two orders of magnitude. This suggests that we can look for a solution which is a double power series in a and k. 
We begin by expanding the solution in powers of a. To first order, we get 



T(t) « 



1 2k, , / c 
(l + 6<5) - 



2k 



(14) 



with 



Tv(r) = 



2k 
37 



(15) 



As seen from the left panel ofFigure4, for k = l/(47r) it isagood approximation in the interval 0.08 fm/c < r < 4fm/c. 
Next, we expand (14 1 in powers of k. In this case, neglecting higher order terms gives 



(16) 



Again, as seen from the right panel of Figure 4, for k — l/47r it is a good approximation whose error is less then 0.9% 
for the temperature evolution. 

IV. OTHER ANALYTICAL SOLUTIONS 

Now we want to show that, with different choices of v, more solutions can be found analytically. These seems to 
be mainly of academic interest. 

Let us first consider the MIT bag model with 



3/3 
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(17) 



Here /? is a positive constant. One can think of this as a consequence of the fact that the model has rj = and ( ^ T. 
The solution to the evolution equation is easily found to be 



Tir) = To 
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1. Note also that (18) decreases with t more 



In fact, it can be added to a list of those in [19] as a special case n 
slowly than (|4]). 

A simple analysis shows that for r < Toe^~'^°'^o /l^ the function (18) grows with increasing t. It gives another 



example of reheating. A way to avoid this problem is to pick tq so that Tq > Toe^~'^'''^o . Together with tqTq = 1/3 
[TB] . this results in the following constraint Tq > 3/3. 
Now let us explore what happens when the fuzzy bag model is accompanied by ^ 



where 7 is a positive constant. In this case, the evolution equation is integrated to give the entropy density 
function of the proper time 



(19) 



as a 



s(r) = 4a-e ^ 
r 



(20) 



^ Curiously, a linear dependence of rj/s on T was derived in |21) within the hard wall model of AdS/QCD. 



7 



where c = tqTq{Tq — ^Q!)e^o . As before, we can obtain an explicit formula for T(r) by solving a cubic equation. It is 



now 
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if re? > J^c 



(21) 



Some comments about formula (21) are in order. Let us first notice that for sufficiently large r this solution 



approaches that of perfect fluid, with c replaced by c. Thus, viscous effects become negligible on long time scales. 
Next, there is apparently a problem of small r. In this case, one must take tq > 7 in order to avoid any unphysical 
reheating. Finally, if we set a = 0, that means a bag model type equation of state, then the solution becomes 



(22) 



The exponential factor on the right hand side shows that the temperature decreases with the proper time more slowly 
in ([22]) than in Q. 



V. CONCLUDING COMMENTS 

There is a large number of open problems associated with the circle of ideas explored in this paper. In this section 
we list a few. 

(i) There is a conjecture that there may be a bound from below on the ratio of the bulk viscosity to the shear 
viscosity [53] 



< 



(23) 



Here Cg is the speed of sound. For the fuzzy bag it is simply = (2T^ — a)/{6T^ — a) yTlj. Combining this with 
(10) and (11), we learn 



T — ka 3m 
6T2 - a - 4^ 



(24) 



Thus, if we accept the idea that the a-corrections are small, we can immediately deduce from ([8|) and (24) that for 
m sa 0.024 the (temperature-independent) ratio 77/s is limited to 



0.08 < 77/s < 0.11 



(25) 



This is a very narrow interval even in the temperature range of interest, where ([ij is valid. Is it reasonable? If not, 
what will happen with the bounds? 

(ii) Given the solution (|6|, we can describe how it behaves for large proper times. ^ To leading order, we have 



T(r) 



1 c 



+ +0(t-^). 



(26) 



This result shows that a small a expansion is no more appropriate. 

Interestingly, the temperature approaches the value \/a/2 which is, according to [5^, the critical temperature.^ 
In other words, there is a non-zero limiting value of T(t) which turns out to be Tc- It is worth mentioning that the 
critical slowing down was also observed in the MIT bag model |19j under the assumption that the bulk viscosity has 



* Clearly, this is of academic interest only. 

® Note that for a = 0.09 GeV^ a simple estimate of a/2 yields Tc = 0.21 GeV that looks quite satisfactory. 
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a power singularity at T = Tc |27j . It seems that this is not the case for regular viscosity coefficients, as seen from 
the examples (|9|, ([l8| and ([2|]). 

The above conclusion requires a caveat. The fuzzy bag pressure ([I]) is negative at Tc- The reason for this is very 
simple and is contained in the formula (jlj. It is valid for T > T„i„, where T,„i„ is close but not exactly equal to T^. 
So, this doesn't mean that cavitation occurs. It means that the model for the pressure (together with those for the 
viscosities) needs to be refined. It would be interesting to see what happens if one does so.^'^ 

(iii) While the lattice data for pure gauge theories suggest that at high temperatures the equation of state of the 
original bag model does get modified by the term [Sl[5], the situation with the equation of state in 2 + 1 flavor 
QCD is far from perfect. There is a discrepancy between the data recently reported by the Wuppertal-Budapest 
collaboration [55] and those of the HOT-QCD collaboration [3]. While this issue is being or will be resolved by 
further lattice simulations, it is tempting to consider a general polynomial model for the pressure 



p(T) = a{T^ + 87T^ - qT^ + 4/3T) - B . (27) 

This is a 3-parameter deformation of the original bag model which reduces to the fuzzy bag ([T]) for /3 = 7 = 0. 

In the absence of dissipative effects, we can obtain an explicit formula for T(r) by solving a cubic equation. For 
small T it is given by 



T(t) = -27 + /3 




(a + 2472)3 



54/2 



(a + 2472 
54/^ 



(28) 



where / = ct~^ — 167"^ — 07 — /? and c = to{Tq + 6"/Tq — \dT^ + /3). The solution ([6]) emerges in the limit /3, 7 0. 

In Figure 5 we give an example of the proper time evolution of the temperature as it follows from ( 28 ) . Wc fit the 
lattice data for the interaction measure [25] to the polynomial model (27 1, as shown in the left panel. This yields 
a = 2.8, 7 = 6.5 MeV, a = 0.09 GeV^ ^ = (0.14 GeV)^, and B = (0.2 GeV)"*." As before, we take Tq = 0.8 GeV and 
To = 0.08 fm/c. What we see from the right panel is that the temperature decreases with the proper time more slowly 
in the presence of the additional terms in the equation of states. Moreover, a difference between the deformations of 
the equation of state becomes visible on the temperature curves when t > 1 fni/c. 



0.5 r 



0.4 
0.3 
0.2 
0.1 



T (GeV) 



0.2 0.3 0.4 0.5 0.6 0.7 0.. 



T (GeV) 



, 0.8 




\ 0.7 




\ 




0.4 




0.3 





0.5 



1.0 



1.5 



2.0 



T (fin/c) 



FIG. 5: Left: The interaction measure normalized by T^. The curve corresponds to (27 1. The dots come from the data [29] at 
R = 28.15 (physical quarfc masses). Riglit: Evolution of the temperature. As in Figure 1, the lower blue and upper red curves 
correspond to the MIT and fuzzy bag models. The middle curve represents ( 28 1 . 



(iv) Here we used the simplified model for a boost-invariant expanding fluid which allowed us to find the solutions 
to the hydrodynamic equations in the first order theory. The main reasons for doing so were to find the solutions 
analytically and to avoid the issue of uncertainty about the transport (relaxation) coefficients at higher orders. 
Although the model is a good approximation to the early time evolution of the QCD matter, it would be very 
interesting to see whether our findings can stand the test of a 3 + 1-dimensional hydrodynamic code. 



A step in this direction was taken recently in 1281 . where the lattice result of [9] was used at low (T < 0.3 GeV) temperatures. 
In doing so, we take the same values of a and B as in the fuzzy bag model. 
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APPENDIX 

The purpose of this appendix is to describe the evolution of the energy density, pressure and trace anomaly in the 
absence of dissipative effects. The corresponding calculations are trivial but a little tedious. 
Combining the solution ^ with the equation of state (II]) , we obtain 




It is instructive to look at the expansion of £(t) in powers of r. To first order, it is given by 



4 2 

£(t) = 3a(^)' +aa(^)' +0(1). (A.3) 

In this formula the second term is the leading a-correction to the evolution of the energy density of the MIT bag 
model [12 



e{r)=3a[^y . (A.4) 

Because it is positive, the energy density of the fuzzy bag falls faster. This explains the curves in the right panel of 
Figure 1. Moreover, we can get an idea about the evolution of the trace anomaly (interaction measure) in this time 
interval. It is 



e — 'Sp 
2^4 



2 

(r) = 2aa(^)' +0(1). 



(A.5) 



Actually, the above analysis can be easily extended to the solution (21 1 by replacing c with ce 
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